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It has been postulated that partonic orbital angular momentum can lead to a significant double- 
helicity dependence in the net transverse momentum of Drell-Yan dileptons produced in longitudi- 
nally polarized p+p collisions. Analogous effects are also expected for dijet production. If confirmed 
by experiment, this hypothesis, which is based on semi-classical arguments, could lead to a new ap- 
proach for studying the contributions of orbital angular momentum to the proton spin. We report 
the first measurement of the double-helicity dependence of the dijet transverse momentum in longi- 
tudinally polarized p + p collisions at y/s = 200 GeV from data taken by the PHENIX experiment 
in 2005 and 2006. The analysis deduces the transverse momentum of the dijet from the widths of 
the near- and far-side peaks in the azimuthal correlation of the dihadrons. When averaged over the 
transverse momentum of the triggered particle, the difference of the root-mean-square of the dijet 
transverse momentum between like- and unlike-helicity collisions is found to be — 37±88 stat ±14 syst 
MeV/c. 

PACS numbers: 13.75.Cs, 14.20.Dh, 21.10.Hw 



I. INTRODUCTION 

Since the startling 1989 result of the European Muon 
Collaboration, which revealed that much less of the pro- 
ton spin is carried by the quark and antiquark spins than 
previously expected [l[, there has been great interest in 
the angular momentum structure of the nucleon. Sub- 
sequent deep-inelastic scattering experiments have con- 
firmed that only ^20-30% of the proton spin is due to 
quark and antiquark polarization [2, 0] . 

The remainder of the spin of the proton must be due 
to gluon spin and/or partonic orbital angular momentum 
(OAM). It is known that the proton anomalous mag- 
netic moment requires orbital angular momentum of the 
quarks, although the combined orbital angular momen- 
tum of all flavors may be close to zero (see e.g. [4, 5]). 
Recent measurements of AG, the gluon spin contribution 
to the proton, are still statistically limited but have ex- 
cluded large values of gluon polarization @, 0, Q , and the 
most recent global study indicates nearly vanishing gluon 
polarization in the presently accessible x range, together 
with a small quark polarization 0. Forthcoming data 
from the Relativistic Heavy Ion Collider (RHIC) should 
place tighter constraints on AG and shed new light on the 
spin puzzle. Meanwhile, progress in the quark and gluon 
helicity distributions has served to help fuel the increas- 
ing interest in orbital angular momentum that began in 
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the 1990s. 

It is important to note that while the total spin of 
the proton as h% is well defined, there is no unique way 
to describe the decomposition of the angular momen- 
tum among the interacting partons within a nucleon (see 
e.g. [13] )■ Thus discussions of partonic orbital angu- 
lar momentum in the proton typically involve a num- 
ber of subtleties, despite the relatively intuitive nature 
of the concept. Two decompositions of nucleon angu- 
lar momentum that have become standard are that of 
Jaffe and Manohar [ll| and that of Ji [12]. While at 
present no quantitative method is known to probe ex- 
perimentally the partonic OAM of the Jaffe-Manohar de- 
composition, in Ji's paper he proposes the experimental 
technique of deeply-virtual Compton scattering to ac- 
cess quark OAM via Generalized Parton Distributions 
(GPDs). Several groups have alreadypursued this ex- 
perimentally challenging path [l3|, [TJ, [l5|, [l|| [TtJ ■ Initial 
measurements of hard exclusive leptoproduction of vec- 
tor mesons, another means of accessing GPDs, have also 
been performed [H, El- Within the Ji decomposition, 
results for the OAM of up and down quarks have become 
available from lattice quantum chromodynamics (QCD) 
calculations [20(. These lattice QCD results suggest that 
the orbital angular momentum for u and d quarks sep- 
arately is quite substantial, but that these contributions 
largely cancel in the proton. 

Another approach to studying the transverse mo- 
tion of quarks and gluons within the nucleon is 
through transverse-momentum-dependent parton distri- 
bution functions (TMDs). The first attempt to use a 
TMD to describe the large transverse single-spin asym- 



metries (SSAs) observed in polarized hadronic collisions 
was made by Sivers in a 1990 paper [2l[, and the 
various TMDs contributing to the leading-order polar- 
ized semi-inclusive deep-inelastic scattering (DIS) cross 
section were laid out by Mulders and Tangerman in 
1996 [23] . Progress was made in both experiment and 
theory throughout the decade, but it was onl y a fter some 
key theoretical developments in 2002-03 [231123. |25| that 
an ongoing period of intense theoretical and experimental 
activity regarding TMDs began. It should be noted that 
thus far, no model-independent quantitative relationship 
between TMDs and parton orbital angular momentum 
has been derived [IE E3 , and it is not clear at present if 
the OAM to which TMDs could provide sensitivity would 
fit within cither the Jaffe-Manohar or Ji decomposition 
of nucleon angular momentum. 

While the majority of work related to investigating 
OAM of the partons within the nucleon has taken place 
since the 1990s, an early theoretical discussion of or- 
bital angular momentum inside hadrons was published 
by Chou and Yang in 1976 [HI, describing the "hadronic 
matter current" inside a polarized hadron. After the 
EMC result [l|, Meng et al. [1^] built upon these semi- 
classical ideas and proposed two experiments to access 
rotating constituents in the nucleon, one in semi-inclusive 
deep inelastic scattering of unpolarized leptons on trans- 
versely polarized protons, and the second in the measure- 
ment of the net transverse momentum of Drell-Yan pairs 
in collisions of longitudinally polarized protons. The lat- 
ter lays the theoretical basis for this analysis: if the trans- 
verse momentum of the partons in the initial state is cor- 
related with the (longitudinal) spin direction, then hard 
collisions involving these circulating partons will lead to 
final states with a net transverse momentum pt with 
magnitude dependent upon the relative orientation of the 
spin directions and the impact parameter of the collision, 
as can be seen in Fig. [T] 

For a particular helicity combination, e.g., positive on 
positive, the transverse momenta of the rotating partons 
add for peripheral collisions and give a net transverse 
momentum to the lcpton pair (in the case of Drell-Yan). 
For small-impact-parameter collisions in the like-helicity 
combination, the helicity-correlated transverse momenta 
of the partons mostly cancel. In the other helicity com- 
bination (unlike-sign) , the opposite effect is seen, i.e., pe- 
ripheral collisions give a small net transverse momentum, 
while small-impact-parameter collisions give a larger net 
transverse momentum. 

The correlation of the parton transverse momentum 
with the orbital angular momentum is expected to de- 
pend on the spatial position of the parton in the proton. 
However, experimentally there is currently no technique 
for determining the impact parameter of an inelastic p+p 
collision, and more specifically the spatial location of the 
parton-parton hard scattering within that geometry. De- 
spite this limitation, in [29j . with a rather simple pic- 
ture of the transverse spatial distribution (homogeneous 
sphere) and momentum distribution (rotational momen- 
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like-sign helicities unlike-sign helicities 




like-sign helicities unlike-sign helicities 




FIG. 1: (Color online) Colliding protons are represented by 
overlapping circles, with proton momentum designated by the 
central symbol, and spin direction designated by the clockwise 
or counter-clockwise arrows. A positive correlation between 
parton transverse momentum and proton spin has been as- 
sumed. For a like-sign helicity combination (positive on posi- 
tive, left panels) , the transverse momenta of the rotating par- 
tons add for peripheral collisions (top left) and result in a net 
transverse momentum of the lepton pair (in the case of Drell- 
Yan) and mostly cancel for small-impact-parameter collisions 
(bottom left). In the other helicity combination (unlike-sign), 
the opposite effect is seen, i.e., peripheral collisions lead to 
cancellations of the transverse momentum (top right), while 
small-impact-parameter collisions give a larger net transverse 
momentum (bottom right). 



turn, kef,, independent of position inside the proton), it 
was found that approximately half of the maximum effect 
((?5y)max = when the vector transverse momenta are 
exactly aligned) remains after integrating over the im- 
pact parameter. This result is based on a semi-classical 
model, with the assumption that all interacting partons 
have the same rotational momentum. As in the case of 
TMDs, there is at present no well-defined relationship 
between the partonic OAM to which this method could 
provide sensitivity and either the Jaffe-Manohar or Ji de- 
composition of nucleon angular momentum. However, it 
is interesting to note that unlike effects due to the Sivers 
TMD [23], in the semi-classical model in which the cur- 
rent analysis is framed, the effect discussed below does 
not require an initial- or final-state interaction to gener- 
ate a non-zero effect. 



II. DRELL-YAN VS. JET k T 

Here, we propose to probe the spin-correlated trans- 
verse momentum of partons within longitudinally polar- 
ized protons involved in hard collisions leading to jet-like 
events at the PHENIX experiment at RHIC. However, in 
PHENIX, due to our limited acceptance (in the central 
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region, A0 = |x2 and \rj\ < 0.35 [3(|), we do not recon- 
struct the true jet kinematics to access the jet transverse 
momentum. An alternative method has been developed 
[3ll ] that examines the dihadron azimuthal angle correla- 
tion to extract the average parton transverse momentum, 
■J (fcj.), on a statistical basis for two subsets of the data, 
like-helicity collisions and unlike- helicity collisions, which 
can then be compared as a measure of the helicity depen- 
dence of the net interacting parton transverse motion. 

Since, in contrast with the Drell-Yan experiment pro- 
posed in we deal here with hadronic final states, 
there could in principle be spin-dependent contributions 
to the measured yj {k'j,} which are not related to the ini- 
tial partonic transverse momentum. The measured di- 
hadron transverse momentum, (p^ ut ), is a convolution of 
the measured fragmentation transverse momentum (jj,) 
and the extracted partons' transverse momenta (k^). 

With the factorization ansatz for the mean p|, of the 
scattered partonic pair presented in (3lj |. 

= = <b e , (i) 

where the superscripts /, S and H denote intrinsic, soft 
(one or several soft gluons emitted) and hard (NLO) con- 
tributions respectively, one might attempt to understand 
the helicity dependence of each term. The conclusion 
of [291 ] is that the difference in the intrinsic contribu- 
tion to the mean square parton transverse momentum 
between positive- and negative- helicity protons, A(fc^) / , 
could be non-zero, since, with a net orbital angular mo- 
mentum, there would be a non-zero helicity difference in 
the vector-summed kr of the initial partons. A(k^) H 
could also be non-zero, e.g., given a helicity dependence 
of three-jet events. This contribution is theoretically cal- 
culable in perturbative QCD, and experimentally, con- 
tributions from a hard component should be accessible 
by measuring and comparing the spin-dependent kx dif- 
ference for several center-of-mass energies. As in QED 
[12], soft radiation in QCD is independent of the polar- 
ization of the emitting particle, so the (fcy) term would 
not contribute to any spin-dependent (kj,) difference. 

Additionally, since (jj) is used to extract (kj,) from 
(Pout); ^ i s important to note that any possible spin de- 
pendence of (jj) can be measured directly in this analy- 
sis. 

The relationship of a measured yj (kj) difference to a 
partonic orbital angular momentum is non-trivial. One 
can attempt to relate the spin-correlated parton trans- 
verse momentum to this difference: 

a<4) j = J2° iiwij {v% • 4> ++ - (% • ( 2 ) 

where the sums are over all partons in the colliding pro- 
tons, is the probability of an interaction of the i th 
and j th partons leading to the final state, is the 



(unknown) impact parameter weighting for the interac- 
tion, and k l T and k 3 T are the two-dimensional partonic 
transverse momenta. In the case with no spin-dependent 
transverse momentum (no orbital angular momentum), 
the difference between the ++ (like-helicity) and the 
H — (unlike- helicity) terms vanishes. The c y can be 
calculated from parton distribution functions, whereas 
the W lJ may be estimated from simulations, given a 
model for the impact-parameter-dependent parton dis- 
tributions. 

It is evident from Eq. [5] that the mixture of initial- 
state partons leading to a 7r° — h final state will have an 
impact on the interpretation of the data. In the central 
arms of PHENIX, where tt° — correlations are mea- 
sured, PYTHIA [33| simulations show that ~ 50% of the 
events leading to ir° — h events are g — g in the initial 
state at ir° transverse momenta below 4 GeV/c (from 4-7 
GeV/c the fraction is ~ 40% ), g— q initial states mak- 
ing the majority of the remainder. Only a small fraction 
of the events are like-flavored q — q in the initial state. 

It is instructive to examine what happens if the sign 
of the orbital angular momentum is different for different 
flavors. When two partons with the same sign OAM in- 
teract in a peripheral p + p collision, then the transverse 
momentum adds constructively as in the top left panel 
of Fig. [H regardless of the sign of the OAM. However, 
if the two interacting partons have opposite sign OAM, 
then the result would be as in the right side of Fig. Q] 
Therefore, an equal mixture of parton interactions with 
like-sign OAM with unlike-sign OAM would result in a 
zero proton-helicity difference in the RMS transverse mo- 
menta. On the other hand, if partons of a certain flavor 
carry no OAM, then interactions involving that flavor 
would contribute nothing to the effect in either helicity 
case, and only act as a dilution to the overall transverse 
momentum difference. 

Given the dominance of gluon scattering for the kine- 
matics of this measurement, then, the results could be 
qualitatively interpreted (within the semi-classical model 
presented) as due to a diluted contribution of the gluon 
orbital angular momentum to the partonic kx- A more 
quantitative interpretation would require a model for the 
OAM dependence on flavor and kinematics, together with 
a process and experimental simulation. 



III. k T FROM DIHADRON AZIMUTHAL 
CORRELATIONS 

In this analysis we used PHENIX high-p^ photon trig- 
gered data from RHIC runs in 2005 and 2006 at y/s = 200 
GeV, as has previously been published for the PHENIX 
7T° cross section asymmetry (All) analysis @, [HI with 
integrated luminosities of 2.5 pb _1 and 6.5 pb _1 respec- 
tively. Neutral pions were selected from photon pairs 
falling in the invariant mass region within M T o ± 2. Oct. 
The signal-to-background ratio for the 7r°'s in the range 
p T = pTt > 3.0 GeV/c is above 15. 
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FIG. 2: (Color online) Azimuthal distributions for real 
(solid curve): dN rea i / dA<f> and mixed event (dashed curve): 
dNmix/ dA(j) pairs. 

The azimuthal correlation function is obtained by mea- 
suring the distribution of the azimuthal (around the 
beam axis) angle difference, A<p = <fi t — 4> a , between a 
7T° (triggered particle) and a charged hadron (associated 
particle). The data is analyzed in eight bins of ir° trans- 
verse momentum from 2.0 GeV/c < prt < 10-0 GeV/c, 



and the associated charged hadron transverse momen- 
tum p\ = pTa bin is selected to be within 2.0 GeV/c < 
PTa < 5.0 GeV/c throughout this analysis. Whenever a 
7T° is found in the event, the real (dN rea i / dAtfi) and mixed 
(dN m - lx /dA(j)) distributions are accumulated. The mixed 
event distribution is applied as a correction factor to ac- 
count for the limited PHENIX acceptance. Mixed events 
are obtained by pairing a ir° taken from a dihadron event 
with many charged hadrons taken from different events, 
randomly selected from a minimum bias data set (no 
high-px photon required) without regard to helicity. The 
mixed event distribution is kept the same for both helicity 
combinations. Figure [2] shows the real and mixed event 
distributions for 3.5 GeV/c < prti^ ) < 4.5 GeV/c. 



The fragmentation transverse momentum \f(Jj) and 
the partonic transverse momentum yj (kj,) are related to 
the widths of the two peaks in the correlation function; 
around A(f> = degrees to obtain a near , and around 
Acf) = 180 degrees to obtain yj (p*|ut) ( tue RMS trans- 
verse momentum of the charged hadrons with respect to 
the 7r°'s). The raw dN rea i/dA(j) distribution is fit with 
the following function to obtain a near (based on a near- 
side Gaussian) and yj (Pout) (based on a more compli- 
cated away-side functional form, as derived in [3ll]): 



dN, 



real 



1 dN miy 



dA4> N dA(j> 



C' + Ci ■ GaUS(0, a ne ar) + C 2 



dN 



far 



dA(j) 



3ir/2N 
tt/2 , 



(3) 



where 



dAcj) 



3-7I-/2 



-p Ta COS A(f) 



,/2 y/27r(pl ut )Erf (V2pr,/V^> 



exp 



Pr a sin2 
2{pLt> 



(4) 



To calculate y/(jf) and y/Jkf) from the cr nmr and 
yj (Pout) values obtained from the fit, the following for- 
mulae from [3ll are used: 



^> = V2. 



PTt ■ PTa 



PTa 



x h 



1 

%h 



r-\/(Pout) 



(il 3/ )(i + 4) 



(5) 
(6) 



where x n = PTa/ PTt, is the analogous ratio of the 
partonic transverse momenta, (z t ) is the ratio of hadronic 
to partonic transverse momentum for the trigger 7r°, and 



(j T y) = » 



Figure[3]and Table[T]show the derived values of (zt) and 
Xh, which were determined through an iterative process 



using a combined analysis of the measured ir° inclusive 
and associated spectra using jet fra gme ntation functions 
from LEP e + e~ measurements [3a, Sal, as in 311. The 



central values were calculated assuming an equal fraction 
of quark and gluon jets, while the systematic uncertain- 
ties on (z t ) and x^ are estimated by taking the RMS 
spread of the g — g, q — q and equal fraction initial-state 
calculations. 



IV. RESULTS 

Fits of the dN rea i / dA<f) distributions were done in three 
ways: 1) all data taken together (summed over spin direc- 
tion), 2) data separated into events from like- helicity and 
unlike-helicity collisions, and 3) the data set randomly 
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FIG. 3: (Color online) Values of derived (z t ) and Xh as ex- 
plained in the text. 



TABLE I: Calculated values of Xhand (zt)for the combined 
2006 and 2005 data sets. 



PTt 

GoV/c 


&h 


M 


2.0-2.5 


1.061 ± 0.003 


0.42 ± 0.05 


2.5-3.0 


0.994 ± 0.000 


0.46 ± 0.06 


3.0-3.5 


0.952 ± 0.004 


0.50 ± 0.06 


3.5-4.2 


0.926 ± 0.008 


0.53 ± 0.06 


4.2-5.2 


0.905 ± 0.011 


0.56 ± 0.07 


5.2-6.5 


0.890 ± 0.014 


0.60 ± 0.07 


6.5-8.0 


0.884 ± 0.013 


0.63 ± 0.07 


8.0-10.0 


0.866 ± 0.013 


0.66 ± 0.06 



separated into two sets of approximately equal number of 
events with the likc-helicity and unlikc-hclicity collision 
type assigned randomly. The first is done as an update 
to our previously published results from 2003 data [3l| 
with higher statistics, but with a slightly different asso- 
ciated charged hadron transverse momentum range, and 
to set the baseline for the partonic transverse momen- 
tum. The second is the measurement of interest, i.e., the 
difference in the net two-parton transverse momentum in 
like- versus unlike-helicity collisions. The results of this 
measurement are to be compared to the model result of 
[29| . The final fitting of the randomly assigned helic- 
ity combinations is done as a measure of the statistical 
accuracy of the fitting results, as explained below. 



TABLE II: Fit parameters a near and yj (Pout) extracted from 
the helicity-averaged 2006 data set. The 2005 results are con- 
sistent within uncertainties. The uncertainties on the param- 
eters do not scale directly with overall statistics, as discussed 
in the text. 



PTt 


(PTt) 


(PTa) 


&near 




GcV/c 


GeV/c 


GeV/c 




GeV/c 


2.0-2.5 


2.23 


2.65 


0.240 ± 0.001 


1.53 ± 0.02 


2.5-3.0 


2.73 


2.67 


0.226 ±0.001 


1.42 ± 0.01 


3.0-3.5 


3.22 


2.71 


0.213 ± 0.001 


1.38 ± 0.02 


3.5-4.2 


3.80 


2.75 


0.199 ± 0.001 


1.28 ± 0.02 


4.2-5.2 


4.61 


2.80 


0.187 ± 0.001 


1.18 ± 0.02 


5.2-6.5 


5.70 


2.86 


0.174 ± 0.002 


1.09 ± 0.02 


6.5-8.0 


7.04 


2.92 


0.167 ± 0.003 


1.00 ± 0.02 


8.0-10.0 


8.78 


2.94 


0.158 ± 0.004 


0.96 ± 0.03 



B. Helicity- Sorted ^(jf^and y/ {k%) 



A. Helicity- Averaged yf (fcy) and y/ (jj,) 

The helicity-averaged fit results are enumerated in Ta- 
ble [ED for the 2006 data set. The results from the 2005 
data are almost identical, with somewhat larger errors. 
The uncertainties on the fit parameters do not scale with 
statistics across the transverse momentum bins, as the 
uncertainty on the extraction of the width of a Gaussian 
distribution which is superimposed on a constant back- 
ground does not scale with the statistics alone, but also 
depends upon the width of the Gaussian. Since the width 
of the peaks depends upon the PTt bin, the uncertainties 
do not scale with the statistics in each bin. Final statis- 
tical uncertainties on the fit parameters are determined 
by a statistical technique discussed below. 

The helicity-averaged y/Jjj) and yj Wr) results for the 
combined running periods are shown in Table Hill and in 
Fig. fH where they are also compared to the previous 
results [3lJ - Note in Fig. [4] that the associated charged 
hadron transverse momentum bin is somewhat higher in 
the current analysis, but when checked by lowering the 
lower limit on pt&, the two results are consistent. 



The process of extracting \/\Jt) and yj (k\) was re- 
peated using the two subsets of the data corresponding 
to collisions involving like- and unlike-helicity protons at 
the PHENIX collision area. Since any spin-dependent ef- 
fects should scale with the polarization of each beam, all 
helicity differences are scaled by p ^ Py , where Pb and Py 
are the run-averaged beam polarizations for the two col- 
liding beams, "blue" and "yellow" respectively, and are 
P B = 0.50 and P Y = 0.49 in 2005, and P B = 0.56 and P Y 
= 0.57 in 2006. Uncertainties on the polarizations were 
propagated together for the two data sets, resulting in a 
4.8% scale uncertainty in the spin-dependent differences. 

The helicity-dependent differences for yj (fcj) and 
yj (jy) (averaged over the 2005 and 2006 data sets) are 
shown in Fig. [5] No yj (j^) difference is observed in any 
PTt bin, and if we assume no pxt dependence and take 
the average over the pxt bins, then the average value of 
the difference in the fragmentation transverse momentum 
is Ay/{j%) = -3 ± 8 stat ± 5 syst MeV/c, consistent with 
zero. 

As discussed earlier, there is no quantitative expecta- 
tion in the difference in Ay/ (k^), but any non-zero mea- 
surement can be attributed to a convolution of initial 
and hard scattering effects. Since no prtdependence is 
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TABLE III: Combined (2005 and 2006) results for U (jf,), \J (fcy) and the helicity-sorted differences. First errors are statistical, 
second are systematic. Statistical and systematic errors are determined as described in the text. 



PTt 


# 








A^/{k 2 T ) 


KjjC V / C 


p(l VPS 


VjC v / c 


VjC v / c 


VjC v / c 


VjC V / c 


2.0-2.5 


792579 


0.582 ± 0.001 ± 0.001 


2.96 ± 0.03 ± 0.36 


-0.008 ± 0.011 ± 0.015 


-0.22 ± 0.19 ± 0.05 


2.5-3.0 


479497 


0.613 ± 0.002 ± 0.001 


2.83 ± 0.03 ± 0.40 


-0.009 ± 0.013 ± 0.015 


-0.11 ± 0.19 ± 0.03 


3.0-3.5 


263174 


0.624 ± 0.002 ± 0.002 


2.87 ± 0.03 ± 0.38 


0.007 ± 0.016 ± 0.015 


-0.17 ± 0.22 ± 0.04 


3.5-4.2 


180554 


0.626 ± 0.003 ± 0.004 


2.79 ± 0.03 ± 0.36 


0.000 ± 0.019 ± 0.015 


0.32 ± 0.22 ± 0.05 


4.2-5.2 


101313 


0.630 ± 0.003 ± 0.006 


2.80 ± 0.04 ± 0.35 


-0.014 ± 0.023 ± 0.015 


0.22 ± 0.24 ± 0.04 


5.2-6.7 


41827 


0.634 ± 0.005 ± 0.009 


2.91 ± 0.05 ± 0.34 


-0.005 ± 0.034 ± 0.015 


-0.03 ± 0.33 ± 0.02 


6.7-8.0 


17916 


0.639 ± 0.008 ± 0.005 


3.01 ± 0.07 ± 0.39 


0.049 ± 0.053 ± 0.015 


-0.36 ± 0.45 ± 0.04 


8.0-10.0 


6775 


0.634 ±0.012 ±0.004 


3.18 ±0.11 ±0.31 


-0.024 ± 0.082 ± 0.015 


-0.48 ± 0.75 ± 0.05 



expected in the model, and the data are consistent with 
a flat distribution, the difference is averaged over prtto 
get A^p^y = -37 ± 88 stat ± 14 s y st MeV/c, consistent 
with zero. 



C. Discussion of Uncertainties 

To check for possible systematic errors due to spin- 
related beam properties or efficiencies, the beam polar- 
ization signs were randomly chosen for each event with 
an equal probability, the A</> distributions were obtained 
for the two false helicity combinations, and the fit param- 
eters extracted. This process was repeated many times, 
giving distributions of the fit parameters that were well 
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FIG. 4: (Color online) Helicity- averaged */ (j^,) and \J (fc|,)for 
combined 2005 and 2006 data. The systematic uncertainties 
on the v/(4>"( green bars) are due mainly to the systematic 
uncertainties on the (zt) and Xh extractions discussed in the 
text. 
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FIG. 5: (Color online) Difference in -J (j'f,) (top panel) and 

yj (fc|,) (bottom panel) for like- minus unlike-helicity com- 
binations. A systematic uncertainty of 4.8% on the verti- 
cal scale due to uncertainty in the beam polarizations is not 
shown. However, this uncertainty only affects the relative 
vertical scale. 



fitted with normal distributions. The widths of the fit- 
parameter distributions for the two false-helicity com- 
binations are then related to the statistical fluctuations 
of the fit parameters. Comparison of these widths with 
the errors returned from the fit indicated that the errors 
on the fit parameters were too small by a maximum of 
~ 15%, especially for p out at the larger pxt bins. 

In order to investigate this non-statistical nature of 
the fit parameter errors, a Monte Carlo simulation was 
employed. Randomly created distributions based on the 
shapes of the real data azimuthal distributions as a func- 
tion of PTt were fitted, extracting the fit parameters and 
errors. This could then be repeated many times, after 
which the widths of the normal distributions of the ex- 
tracted fit parameters were compared to the fit errors. 
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The exact same trend as a function of prt was seen in 
the Monte Carlo - the fit parameter errors were underes- 
timated at the larger values of prt by ~ 15%. Since the 
Monte Carlo is purely statistical, these results reflect the 
true measure of the statistical uncertainty of the fit pa- 
rameters. The statistical uncertainties presented for the 
data are thus those obtained from the spin-randomization 
procedure described above. 

The dominant systematic uncertainties were deter- 
mined from the uncertainties in £h and (zt). Addition- 
ally, the distribution of fit parameter values for different 
mixed-event distributions was fitted with a normal distri- 
bution, the width of which was a measure of the system- 
atic uncertainty of the fit parameters due to the detector 
acceptance correction process. 



V. DISCUSSION 



The smallness of A -J (jj) confirms the expectation 
that transverse momentum effects in the fragmentation 
should not be large in processes with a longitudinally po- 
larized initial state and thus simplifies the interpretation 

of A v /(fcfy. 

Comparing our measured value of Ay/ (fc|,) to the cal- 
culation of [291 ] . 

A(p 2 T ) » 1.9<fc ) 2 (7) 

yields: 

« A\f^) = -37 ± 88 stat ± 14 syst MeV/c, (8) 

approximately an order of magnitude less than the par- 
ton intrinsic transverse momentum associated with the 
uncertainty limit. Assuming all contributions to this dif- 
ference come from intrinsic parton motion, and taken to- 
gether with the expected level of contribution from the 
g~g channel and our model assumptions, this could qual- 
itatively suggest a small gluon orbital angular momentum 
in a longitudinally polarized proton, integrated over our 
kinematic region. A more direct connection between this 
measurement and partonic OAM is complicated by the 
subprocess contributions and unknown impact parameter 
and transverse position space weighting of the partons. 
In addition, further theoretical work is needed to place 
the model of [2!| within a rigorous QCD framework. We 
hope that the measurement presented here will serve to 
encourage the theory community to pursue this task. 

As discussed in Section HI since the 1990s there has 
been intense interest in partonic OAM, or more gener- 
ally, in the non-collinear motion of partons within the 
nucleon, and there are several approaches currently be- 
ing used to attempt to increase our understanding of the 
role that this partonic motion plays in nucleon structure, 
not all of which can be directly related to one another. 



The measurement presented here, inspired by the pro- 
posal in [29j for Drell-Yan production in longitudinally 
polarized p+p collisions, but utilizing a dihadron correla- 
tion technique, represents a novel experimental approach 
to probing partonic OAM. 

A dijet correlation technique in single-transversely po- 
larized p + p collisions has already been used at RHIC to 
probe the Sivers TMD [37| . following a proposal in [38|. 
The current measurement has the potential to probe par- 
tonic OAM in a longitudinally rather than transversely 
polarized proton. Dijet and dihadron correlation mea- 
surements in (polarized) p + p collisions provide an im- 
portant tool to investigate the non-collinear motion of 
partons within the (polarized) nucleon, and ideas for ex- 
panding on the existing techniques would be most wel- 
come. 



VI. SUMMARY 



In conclusion, yf (jj,) and y (fcf.) have been extracted 
from dihadron azimuthal angular correlations in longi- 
tudinally polarized p + p collisions at y/s — 200 GeV. 
The helicity differences for both quantities are consistent 
with zero when averaged over the ir° transverse momen- 
tum range accessible, with a magnitude less than 5% of 
the corresponding spin-averaged quantities. Comparison 
to a similar measurement that can be performed on longi- 
tudinally polarized p+p collisions at y/s— 500 GeV is ex- 
pected to provide additional information regarding hard 
vs. intrinsic contributions to the measured A^/ (fcj.). The 
PHENIX collaboration collected such a data set in early 
2009. Future data at RHIC will increase the statistical 
significance, and upcoming PHENIX upgrades will allow 
measurements in different kinematic regimes, changing 
the partonic mix probed. In the longer-term future, the 
accumulation of large luminosities for polarized p+p col- 
lisions at RHIC should also make possible a Drell-Yan 
measurement, as originally proposed in 29]. 
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